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ABSTRACT 

Expanding on an earlier work that relied on linear force-free magnetic fields, 
we self-consistently derive the instantaneous free magnetic energy and relative 
magnetic helicity budgets of an unknown three-dimensional nonlinear force-free 
magnetic structure extending above a single, known lower-boundary magnetic 
field vector. The proposed method does not rely on the detailed knowledge of 
the three-dimensional field configuration but is general enough to employ only a 
magnetic connectivity matrix on the lower boundary The calculation yields a 
minimum free magnetic energy and a relative magnetic helicity consistent with 
this free magnetic energy. The method is directly applicable to photospheric or 
chromospheric vector magnetograms of solar active regions. Upon validation, it 
basically reproduces magnetic energies and helicities obtained by well-known, but 
computationally more intensive and non-unique, methods relying on the extrap- 
olated three-dimensional magnetic field vector. We apply the method to three 
active regions, calculating the photospheric connectivity matrices by means of 
simulated annealing, rather than a model-dependent nonlinear force-free extrap- 
olation. For two of these regions we correct for the inherent linear force-free 
overestimation in free energy and relative helicity that is larger for larger, more 
eruptive, active regions. In the third studied region, our calculation can lead to 
a physical interpretation of observed eruptive manifestations. We conclude that 
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the proposed method, including the proposed inference of the magnetic connec- 
tivity matrix, is practical enough to contribute to a physical interpretation of the 
dynamical evolution of solar active regions. 

Subject headings: Sun: atmosphere — Sun: corona — Sun: coronal mass ejec- 
tions — Sun: flares — Sun: magnetic fields — Sun: photosphere 



Introduction 



Several decades have elapsed since the notion of magnetic helicity and its application 
to solar magnetic fields were introduced. This considerable time has been marked by an 
impressive volume of published works on the subject. Yet, we are still lagging behind in 
understanding, first, how to practically calculate magnetic helicity in the Sun and, second, 
what is the actual role of magnetic helicity in solar eruptions. 

Undoubtedly, the apparent lack of a breakthrough in this topic stems from our cru- 
cially incomplete knowledge of solar magnetic fields: unable to observe their generation at 
the solar interior and to precisely measure them in the solar atmosphere, we can routinely 
detect and measure them only in the photospheric and/or low-chromospheric interface. To 
calculate magnetic helicity, however, we need either the three-dimensi onal magnetic field in 
all or part of the solar coronal volume, inc luding; its lower boundary (jWoltjerl Il958l : iBerger 



984: 



Finn fc Antonsen 



1985 



Bergerlll999r) or the flow velo city field on this lower boundary 



(IBerger &: Fieldll 19841 : IChaell200ll ; iDemoulin fc Bergerll2003l ) . Both defy a unique calculation 
showing critical model-dependent ambiguities eithe r when extrapolating f or th e coronal mag- 



netic field (see the comprehensive comparisons of ISchrijver et al. 



( 200 81)) or when inferring a reliable photospheric velocity fie 



Kusano et al 



2002 



iNindos &: Zhang 



2002 



Nindos et al 



2003 



hood) and 



Metcalf et al. 



d ( November fc Simo 



Longcopel 12 004 



Schuck 



1988; 



2005. 



2008 



Georgoulis fc LaBontell20()6l : IvVelsch et aDl2007l : iRavindra et aDl2008l : Ichae fc Sakurai 



20081 ). Knowledge of the three-dimensional coronal magnetic field of, say, an active region, is 
needed to calculate the instantaneous magnetic helicity (and energy) budgets in the region, 
while knowledge of the local flow field at the lower atmospheric boundary is necessary to 
calculate the injection rate of magnetic helicity in the solar atmosphere due to the region's 
temporal evolution. The accumulated helicity budgets are then obtained by integrating 



the h elicity injection rate in time (e.g., iGreen et al.ll2002t INindos et al.ll2003l ; lLaBonte et al. 
20071 ). 
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When the role of magnetic helicity in solar eruptions is examined, there are stud- 
i es suggesting that helicity is not necessa ry for flares and coronal mass ejections (CMEs) 
( IPhillips et al.ll2005l ; IZuccarello et al.ll2009[ ). At the same time, other stu dies suggest that he- 
lical ( or helical by proxy) a c tive regions tend to be the most erupti ve ones ( INindos &: Andrews 
20041 ; lLaBonte et aD 120071 : iNindosI bood beorgoulis et al.l l2009h . One should acknowledge 
the fact that solar eruptions can occur, at least in models, in the absence of significant mag- 
netic helicity accumulations. However, since helicity is a signed quantity with right-handed 
(positive) and left-handed (negative) senses, absence of a significant helicity budget could 
also mean significant helicity accumulation of both senses at roughly similar amounts - this 
can le ad even to the so-c alled "helicity annihilation" that is a proposed eruption mecha- 
nism (IKusano et al.ll2003l ). Moreover, several eruption mechanisms stem from instabilities 
that do not explicitly rely on helic ity. Such mecha nisms are th e magnetic flux cancellation 
(Ivan Ballegooijen fc Martenslll989l): "ho op" force rtchenl Il996h : bre akout (Untiochos et all 



19991 ): tether-cutting jMoore et al.lboOll ): and the torus instability kliem fc Toroklbood ). 
among others. On the other hand, a popular eruption mechanism that re l ies on magnetic 
helicity is the helical kin k instability (IRust fc Kumarlll996l ; iBaty et al.lll998l ; iTorok fc Kliem 
20051 ; Kliem et al .1120121 ). Moreover, it has been shown that injection of helicity in a modele d 
eruption results in faster CMEs after a helicity threshold is exceeded ( jZuccarello et al.ll2009l ). 



In theory, therefore, solar eruptions can occur with or without a significant magnetic 
helicity budget, namely without a dominant helicity sense, although eruptions may be nec- 



essary to diffuse into the heliosphere the excess helicity produced in the Sun (iLowl 11994 



Rustl ll994j, l2003l ) . This is because helicity is known to be roughly invariant in a volume 
enclosing an isolated m agnetic struct ure even under resistive manifestations such as mag- 
netic reconnection (e.g.. iBergerl ( 119991 ) and references therein). On the other hand, virtually 
all eruption mechanisms, regardless of helicity dependence, result in strongly helical erup- 
tion ejecta widely known as flux ropes. As to the pre-eruption situation, we cannot collect 
clues about the actual role played by helicity unless we (i) compare the helicity budgets be- 
tween non-eruptive and eruptive active regions in the pre-eruption state, and (ii) assess the 
relevance of proposed eruption initiation mechanisms with observations in a very detailed 
manner. 

A prerequisite of both objectives above is the practical and reliable calculation of mag- 
ne tic helicity in observed s olar m agnetic structures. A first step in this direction was taken 
by lGeorgoulis fc LaBontd (120071 ) - hereafter Paper I. In Paper I we described a methodology 
to simultaneously calculate the relative magnetic helicity and the free magnetic energy, with 
respect to a potential-field reference, of a magnetic structure represented by a single solar 
vector magnetogram. The underlying assumption was the validity of the linear force-free 
(LFF) field approximation in the magnetic structure. Calculation of both the relative helic- 
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ity and the free energy was physically consistent and did not rely on a prescribed flow velocity 
field or the detailed three-dimensional coronal field above the lower-boundary magnetogram. 
In essence, it was a convenient surface calculation of physical quantities that stem from 
the three-dimensional magnetic field on and above the surface. The main drawback of the 
methodology, however, was its central assumption of a constant-a, LFF magnetic structure 
- an as sumption that is known to be unrealistic in both the solar surface and the overlying 
corona (IMetcalf et al.lll995t iGeorgoulis fc LaBontell2004l ; ISocas-Navarroll2005l ). Nonetheless, 
comparison between a non-eruptive and an eruptive solar active region revealed, even beyond 
large uncertainties inherent to the LFF field approximation, that the most profound differ- 
ences between the two regions occurred in their budgets of free energy and relative helicity: 
for a factor of ~ 3-difference in unsigned magnetic flux between the two regions the energies 
and relative helicities were different by a factor of ~ 9, with the largest values assigned to 
the eruptive active region. As we show in this work this very large difference is partly due 
to the adopted LFF field approximation. 

In Paper I we explicitly stated that the proposed methodology would serve as the basis 
for a more realistic, nonlinear force-free (NLFF) field approximation in calculating the mag- 
netic energy and relative magnetic helicity. We take this step in this work. The analysis of 
Paper I is extende d to derive the s e lf term s of free energy and relative helicity, while we draw 
from the study of iDemoulin et al.l (120061 ) to derive the mutual terms of these quantities. By 
construction, the LFF field methodology of Paper I treated a given magnetic structure as a 
single, isolated, force-free flux tube and hence it was unable to predict mutual energy and he- 
licity terms occurring due to the interaction between different flux tubes. This work assumes 
a collection of discrete, slender force-free flux tubes with variable force-free parameters and 
hence calculates both self and mutual terms of energy and helicity. As in Paper I, this NLFF 
field approach is a surface calculation that does not use three-dimensional field extrapola- 
tions or velocity fields. Instead, the proposed method uses a magnetic connectivity matrix 
on the boundary where the vector magnetogram is obtained. This matrix can be obtained 
in any way possible, be it a field extrapolation or not. Therefore, our method is general and 
applies to any connectivity matrix, regardless of inference. To provide perspective, we apply 
the method to the same active- region magnetograms as in Paper I and compare the results. 

The study is structured as follows: the methodology of the calculation is given in Section 
[2j The adopted validation procedure and its results are given in Section [3j Section @] provides 
the numerical results obtained by applying the method to three different solar active regions. 
Section \5\ discusses our findings and provides our conclusion and future perspective. 
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2. Methodology 

2.1. Magnetic connectivity matrix and a-distribution 

The first task is to translate a continuous vector magnetogram into a collection of 
discrete force-free flux tubes with known footpoints, flux contents, and different force-free 
parameters a. If the three-dimensional coronal magnetic field configuration was available, 
then one would be able to trace each magnetic field line separately (here the footpoint of 
a "field line" is restricted to the resolution element [pixel] of the studied magnetogram). 
The coronal configuration may be assessed by extrapolations of various sophistication levels 
(i.e., current-free, LFF- or NLFF-field) but the true configuration is unknown. Moreover, 
tracing and analyzing each field line separately would be impractical and unnecessary. For 
this reason we simplify the vector magnetogram into a collection of thin flux tubes as follows: 



1. We translate the magnetic field configuration into an ensemble of "magnetic charges" 
us ing; the flux partitio n method introduced in the magnetic charge topology model 



of iBarnes et al.l (j2005l ). This is a flux tessellation scheme that relies on a modified 
downhill-gradient minimization algorithm with certain provisions about saddle points. 
This step requires only the normal (vertical) magnetic field component B z . The chosen 
thresholds for partitioning a magnetogram for this work are (i) a threshold of 50 G 
in \B Z \, (ii) a minimum mag netic flux of 10 20 Mx per partition, and (iii) a minimum 
area of 40 magnetogram pixels per partition. These criteria are set to prevent the 
inclusion of quiet-Sun, weak-field, and very small-scale structures, respectively, into 
the calculation, unnecessarily adding to both complexity and required computing time. 
Only partitions that satisfy all three threshold criteria are selected for further analysis. 
Upon completion, we can readily assess the flux content and flux-weighted centroid 
position of each magnetic flux partition. 

2. Assuming that flux partitioning returned p positive-polarity and n negative-polarity 
magnetic partitions, together with their respective fluxes J^; i = {1, ...,p} and J^j; 
j = {1, ...,n}, one may define apxn magnetic-flux connectivity matrix. The matrix 
will contain the fluxes ^f° n committed to the connection ij between the i-positive- 
polarity and the j-negative-polarity partition. Obviously, ^f° n = in case the two 
partitions are not connected. Along with the flux connectivity matrix we construct one 
more p x n matrix containing the vector positions of the two flux-weighted centroids 
of connected partitions. 

3. Each magnetic connection is hereafter assumed a slender flux tube with flux con- 
tent &i° n and footpoints corresponding to the flux-weighted centroids of the two in- 
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volved partitions. To determine the force-free parameter a of this tube we find the 
a-parameters for each partition. From the force-free approximation one may easily 
deduce that the flux-weighted mean a-value over a magnetic partition of flux & is 
given by 

4tt I 

where / is the total electric current of the partition and c is the speed of light. The 
total current I can be calculated by using the integral form of Ampere's law on the 
lower boundary magnetic-field vector B, i.e. 

B-tfl , (2) 

where integration occurs along the closed contour °€ surrounding the partition. 

On the practical side, a valid question is how to determine the bounding contour ^ 
of the partition in order to evaluate Equations (jTJ), (j2J). The partition shapes cannot 
be modeled easily since a partition can assume any closed-curve shape without restric- 
tion. To determine and its contiguous order of points we have developed an "edge 
tracker" that minimizes the length of the curve bounding the partition. Minimization 
is performed by iteratively choosing pairs of neighboring boundary points. This is 
a classic al optimization p roblem that we solve iteratively via a simulated annealing 



method (IPress et al.lll992l ). 

Let Oj, atj be the calculated force- free parameters of the two partitions i and j. We 
assign a force-free parameter 

Oiij = -{oti + atj) , (3) 

for the resulting connection. For each of the two a- values a^, otj there are respective 
uncertainties Sat Saj due to the uncertainties SI in the calculation of the total current 
I (Equation [1]), assuming that the magnetic flux & is known without uncertainty. The 
respective uncertainty Sa^ is, then, 



Sa l3 = ^Saf + Sa] . (4) 

The flux connectivity matrix J£^ on described, we now discuss how we populate it. Ob- 
viously, the result of any magnetic field extrapolation can be translated into a connectivity 
matrix by tracing all extrapolated field lines that open and close within the lower boundary. 
At this stage, we ignore magnetic connections closing beyond the limits of the finite lower 
boundary. Tracing closed field lines from footpoint to footpoint, we add the flux contents 
of field lines that are rooted in the same pair of partitions, thus constructing JPf° n . The 
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simplest connectivi ty matrix, ^/f" t , is the one obtai ned by a current-free (potential) field 
extrapolation (e.g.. ISchmidtl 1 1 964l ; lAlissandrakia 1 198 1[ ) . Any non-potential extrapolation can 
also be used here, but if we use a NLFF field extrapolation we will reach a non-unique result 
subject to the details of the extrapolation metho d. For this reason, o ur m ethod of choice 
is the simulated annealing method introduced by iGeorgoulis fc Rustl (120071 ). The method 
minimizes the magnetic flux imbalance simultaneously with the separation length (footpoint 
distance) of the chosen flux tubes thus emphasizing connections between tightly arranged en- 
sembles of flux partitions, most notably in active regions with pronounced magnet ic polarity 
inversion lines. We have revised the original concept of IGeorgoulis fc Rustl (120071 ) to (i) in- 
clude a mirror flux distribution with as much positive- and negative-polarity magnetic flux as 
the negative- and positive-polarity magnetic flux of the original magnetogram at large (more 
than twice the diagonal length of the original magnetogram) distances, thus producing an 
exactly flux-balanced magnetic structure and treating large-scale, "open" magnetic connec- 
tions, and (ii) include a constant normalization length R m ax equal to the largest length scale 
of the enlarged, flux-balanced magnetogram. These revisions result in a unique connectivity 



matrix ^f° n for the chosen minimization functional 



M 



X) 



R r , 



+ 



(5) 



The above revisions to the original s imulated annea li ng sch eme of IGeorgoulis fc Rust 



( 120071 ) alleviate the criticism applie d bv iBarnes fc Lekal (120 08). First, these authors ar- 
gued that the connectivity result of Georgoulis fc Rust ( 20071 ) was not unique, depending 
on the origin of the coordinate system used, because R ma x in Equation ([5]) was originally 
|rj| + \vj\. Although tests with different system origins showed little, if any, impact for the 
resulting conne c tivity , the introduction of the fixed Rmax puts this issue to rest. Moreover, 
Barnes fc Lekal ( 120081 ) claimed that simulated annealing yields an unphysical connectivity 
matrix that matches neither the potential-field connectivity nor the true coronal connec- 
tivity. Due to our inability to measure the three-dimensional magnetic field vector in the 
corona, however, the "true" connectivity is unknown. Therefore, one cannot comment on 
its similarity, or difference thereof, with the connectivity revealed by simulated annealing. 
We continue to rely on annealing because it emphasizes connectivity in tightly organized 
active regions, that are statistically the most eruptive ones. Point taken, the methodology 
discussed here is more general and can accommodate any connectivity matrix. 

In Figured] we show an example of connectivity calculat ion in NOAA active region (AR) 
10254 , recorded by the Imaging Vector Magnetograph (IVM; lMickey et al.lll996l ; lLaBonte et al 
19991 ) on 2003 January 13. The difference between the potential-field and the simulated- 
annealing connectivity is obvious with the latter committing more flux to fewer, more closely 
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seated, partitions. Also shown is the map of the flux-weighted a-value of each flux partition 
(Equation ([!])). 



2.2. Magnetic Energy and Relative Magnetic Helicity Budgets in the NLFF 

Field Approximation 

Consider a set of N discrete magnetic flux tubes in force-free equilibrium. The magnetic 
helicity of this set can be viewed as the sum of all terms present in a diagonal matrix N x N. 
Diagonal terms I — m (l,m = {1,...,N}) correspond to self -helicity terms and are due to 
the helical features of each flux tube independently. Off-diagonal terms / 7^ m are due to the 
interaction between pairs (/, m) of flux tubes and correspond to mutual-helicity terms. For 
an open volume, w here the set of flux tube s permeates a lower boundary and extends in the 
half space above it, iDemoulin et al.l (120061 ) showed that the relative (with respect to that of 
a potential field) magnetic helicity of the set can be written as 



N 



N 



N 



E^f+E E c. 



(6) 



1=1 



1=1 m=l,l^m 



The two terms of the rhs of Equation ([6]) correspond to the self and mutual helicity of the set 
of flux tubes, respectively. Here T\ is the self-helicity factor of the flux tube I with a magnetic 
flux content $; and is due to its internal structure (twist and writhe). Furthermore, £[ m is 
the mutual-helicity factor due to the interaction of a given pair (I, m) of different flux tubes. 
For the studied open volume, IDemoulin et al.l ( 120061 ) further found 



C, 



lm 



pclose _i_ parch 



(7) 



where Cf^ se is the Gauss linking number, a signed integer reflecting the nu mber and sense of 



the turns a flux tube I winds around a flux tube m and vice versa (see also Moffatt &: Ricca 



(119921 )). and is the mutual-helicity factor of two arch-like flux tubes that are not 



winding around each other. This factor is a real number and can be attributed to the 
translational motions needed to bring the tubes from infinity to their prescri b ed fo otpoint 
positions. Derivation of various £^ c?t - values was described by IDemoulin et al.l (120061 1 and is 
further explained in Appendix A, where additional cases pertinent to our analysis appear. 

Further assuming slender flux tubes, i.e. flux tubes with typical diameter that is much 
smaller than their footpoint separation, to make use of the outcome of the connectivity 
matrix of Section 12.11 the mutual-helicity ter m of Equation (j6[) desc ribes the mutual helicity 
of the set calculated by the flux-tube axes (IDemoulin et al.l 120061 ). However, both T\ and 
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Cim are unknown unless knowledge of the sub-photospheric closures and coronal shapes of 
the flux tubes is available. 

After defining the relative magnetic helicity of t he set of flux tubes, it is necessary 
to define the corresponding free magnetic energy E c . iDemoulin et al.l (120061 ) provided an 
expression for E c only in the case of a closed volume, where the entire length of the closed 
flux t ubes is known and visible. This expression cannot be used here. This being said, Paper 
I and lBergerl (119881 ) defined an energy-helicity formula in the NLFF field approximation. The 
free magnetic energy in this formula reads 



in 



«A • BdV 



(8) 



where V is the integration volume. Obviously AB, where A is the vector potential (V x A = 
B), reflects the volume density of the relative magnetic helicity provided that A obeys the 
Coulomb gauge (Paper I). In the simplified, constant-a (LFF field) approximation, Equation 
(jH) reduces to E c = a/(8ir)H m . From Equation (jSJ), but also from the necessity to ensure 
that E c — >■ when H c — > 0, we approximate the free magnetic energy in our set of N discrete, 
slender flux tubes by the expression 



N N N 

1=1 1=1 m=l,ly£m 



ai£i m $i$ r 



(9) 



thus employing different force-free parameters ol\ for different flux tubes I. Besides ensuring 
self-consistency, Equation (Q makes sure that E c — > faster than a, because in case of a 
potential magnetic field, T) tends to zero and £z m $z<& m -terms algebraically sum to zero, at 
least in a flux-balanced magnetic structure (see Section [2.31) . That both E c and H m must 
tend to zero for potential fields is also a necessity: to view this simply, consider the LFF field 
approximation again where H m ~ E c /a. In case of nearly potential fields, where a — > 0, 
if E c tends to zero slower than a then \H m \ — > oo. We have shown in Paper I that indeed 
H m — > when a — > 0. 

Its advantages given, a we akness o f Equa ti on (ED is that it i s quali tatively similar to the 
current- channel description of iMelrosd (120041 ) . IDemoulin et al.l ( 120061 ) argued convincingly 
that this description is not equivalent to the flux-tube description attempted here. This 
is because each existing flux tube I (/ = {1,...,N}) should spawn a number of additional 
potential flux tubes beyond it in a space-filling, force-free configuration. These additional 
flux tubes induce additional terms in Equation (JHD when interacting with the non-potential 
flux tubes of the set. Otherwise put, for Equation (jUJ) to be valid, the field should not be 
space-filling, that is, "sheath" currents should contain each of the N flux tubes that become 
then embedded in a field-free space occupied by non-magnetized plasma. This might be 
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valid in the photosphe re but is not the case in the corona (e.g., lLongcope fc Welschl 12000 
Georgoulis et al.ll2012l ). where fields are thought to be force-free. Therefore, Equation (Q is 
not fully consistent with the NLFF field approximation that we pursue in this study. 

Despite shortcomings, however, Equation fl9]) can still serve as a lower limit of the 
magnetic free energy. This is because the additional terms induced by the spawned potential 
flux tubes should only add positive increments to E c . Validating our results in Section [3j we 
show that E c from Equation is indeed a realistic lower limit for the magnetic free energy 
of the NLFF field. 

The description of magnetic energy budgets of the set of flux tubes is complete when 
the reference (potential) energy E p is calculated. This can be done in more than one ways, 
by simply using the normal field component B ? at the anchoring , boundary plane of the flux 
tubes. In particular, using the Virial theorem (jMolodenskylll974l ; lAlyill984j ) in its frequently- 
used form (see Paper I for a discussion), we have 



E„ 



— I r • B n B,dS 

47T 



(10) 



where r is a vector position with arbitrary origin in the anchoring plane S and B p is the 
potential-field vector on S. Moreover, by means of the expression 

1 

47T 



E„ = — 



Bp x A p • zdS 



111 



derived in Paper I, where V x A p = B p and z is the upward-pointing unit vector normal 
to the plane S. For a reliable calculation of E p on S one needs a magnetic structure as 
flux-balanced as possible. Finally, by volume-integrating the energy density B^/(8tt) of the 
potential field B p . In Paper I all methods were shown to provide nearly identical E'p-values. 
Then, the total magnetic energy E t of the studied magnetic structure becomes 



Et = E v + E c 



(12) 



In conclusion, the common, self-consistent definition of H m and E c will allow us to derive 
both quantities from observed solar vector magnetograms without requiring sub-photospheric 
or coronal information. This practical way of calculating NLFF field free-energy and relative- 
helicity budgets is described in the following. 



2.2.1. Self terms 



The self terms of the magnetic free energy and the relative magnetic helicity refer 
exclusively to the twist and writhe of each magnetic flux tube independently and are given 



- 11 - 



by 

N N 

E Cs Bl f=e- y E a l T l dOSe ^ and ( 13 ) 
71 1=1 1=1 

respectively. Both terms will be calculated, not by the above formulas, but by generalizing 
the linear analysis of Paper I. 

For the free energy of a single force-free flux tube we derived in Paper I the linearized 
expression 

E c = J r u n d 2 a 2 E p , (14) 

where a is the unique force-free parameter, d is the linear size element (the pixel size, in 
case of observed magnetograms) , and Tu n is a linearized scale factor calculated in Fourier 
space. Assuming a collection of force-free flux tubes rather than a single tube, Equation 
ffl4"j) describes the self term of the free energy under the condition that a is fixed, i.e. the 
same for all tubes. For different a-values per tube, that is, in case of the NLFF field 
approximation, Equation f[T4"j) would need to be evaluated for each tube. This is untenable, 
however, because the potential energy E p (Equations (fTUl) . (TTTT) ) and Tu n are calculated once, 
for the entire plane S of the magnetogram and not for each flux tube separately. To address 
this problem, thus generalizing the analysis of Paper I, we investigate the relation between 
the "scaled" potential energy Ti in E p and the magnetic flux $, in case of a single flux tube, 
or a magnetogram envisioned as a single flux tube. For a large number of magnetograms 
this relation is shown in Figure [2j The straight line represents the least-squares best fit and 
reveals a robust power-law scaling of the form 

T Un E p = A§ lx , (15) 

where A = 1 o~ 16 - 731±0 - 08 and A = 1.153 ± 0.002 are the scaling constants. The substantial 
statistics of Figure [2] stem from the fact that, to calculate J-^ n , E p , and $, one does not 
need vector magnetograms; therefore we have used tens of th ousands of active-reg ion mag- 



netograms recorded by the Michelson Doppler Imager (MDI; IScherrer et al.lll995l ) onboard 
the Solar and Heliospheric Observatory (SoHO) mission. All these magnetograms were lo- 
cated within ±30° from the central solar meridian; hence, the line-of-sight fie l d cor responds 



to the vertical field component within 15%. The correction of iBerger fc Lited (120031 ) was ap 



plied to compensate for the underestimation of the magnetic fluxes in these magnetograms j. 
Substituting Equation ffl5|) into Equation ffT4l) for the case of a single force-free flux tube we 



^Notice that the SoHO/MDI database we possess was constructed before the creation of the Level 1.8 
MDI database. Therefore, underestimation of magnetic fluxes in sunspots and plage still exists in these data. 
Point taken, omitting the correction gives almost identical scaling constants A and A. 
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find 

E c = A{ad)H 2X . (16) 

The key assumption hereafter is that the scaling relation of Equation ffl5|) holds for individual 
flux tubes embedded into a collection of discrete flux tubes with different a-values. Indeed, 
assuming that each magnetogram in Figure [2] corresponds to a flux tube with flux content 
$ (the unsigned magnetic flux of the magnetogram) and force-free parameter a (inferred 
by the linear force-free approximation applied to the magnetogram), Equation ( [16]) would 
provide the free magnetic energy of the flux tube. This energy would be independent from 
the free energy of another flux tube if mutual effects are ignored. Under these conditions, 
we generalize Equation ( I16p for a collection of N flux tubes to provide the self term of the 
free magnetic energy of the ensemble: 

TV 

E Cself =Ad 2 Y,*i$i X • (17) 
i=i 

Equation f lTTI) can replace the respective Equation f fl3|) for E Caa[f , hence accounting for the 
unknown Tf lose in each flux tube I. 

Regarding the self term of the relative magnetic helicity (Equation f lT5]) ). the linear 
force-free approximation of Paper I gives 

H m = Six Fi^d 2 aE p . (18) 

Applying the scaling relation of Equation ffT5]) we find 

H m = 8nd 2 aA$ 2X , (19) 

which is, again, independent for each flux tube if mutual terms are neglected. Assuming 
that this scaling holds for individual flux tubes, we generalize Equation (Tl9|) for a collection 
of flux tubes to obtain the self-helicity of the ensemble: 

N 

H mself = Md 2 A • (20) 

i=i 

This can also replace the respective expression for H mself in Equation ffT3|) . thus accounting 
for the unknown T) in each flux tube /. 

2.2.2. Mutual terms 



By construction, the linear analysis of Paper I cannot be used to calcu late any mutual 
energy or helicity terms. For these terms we will implement the analysis of IDemoulin et al. 
(1200 A . 
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Combining Equations (J7|) and (Q, the mutual term of the magnetic free energy in our 
set of discrete, slender flux tubes is given by 

N N 

Ec mut = ^J2 E atCiT + £%?)*l*rn ■ (21) 
1=1 m=l,lj^m 

As a first step to calculate E Cmut we acknowledge that each term contributed by the inter- 
action of a pair (l,m) of flux tubes provides a positive increment AE Cmut for E Cmut . This 
increment is given by 

AE Cmut = i-(a, + a m )(£^ e + £%?)Q& m , J,m = {1, ...iV} . (22) 

To provide Equation ( 122]) we have also taken into account that £/ m = C m h an d this is the 
case for both Cf^ se and C^ h and for every (l,m). This simply means that the N x N 
helicity matrix is symmetric. Further, we know from theory that Cf^ se = 0, ±1, ±2, ... and 
\Cf^ ch \ < 1 (Appendix A). In our pursuit for a minimum magnetic free energy, therefore, we 
assume that the interacting flux tubes do not wind around each other, so Cf^ se = for every 
(l,m). This leads to a free-energy increment per interaction 

AE Cmut = ^(a l + a m )CZ h ^ m , Z, m = {1, ...N} . (23) 

87T 



According to iDemoulin et al.l (120061 ) and Appendix A there are only two possible values 
of £?m h i n every possible case of interaction for a pair (Z,m) of flux tubes. These values 
depend on whether flux tube I is " above" flux tube m (Cf^ h = £f™ h ) or vice- versa {C^ h = 
C^ h ). The possible cases of flux-tube interaction (see Figure [T]) are (i) intersecting footpoint 
segments of the interacting flux tubes, (ii) non-intersecting footpoint segments (there the two 
values coincide, i.e., = £flh h ), and (iii) "matching" footpoint segments. A discussion of 
"matching" footpoints, where we assume that these footpoints are located within unresolved 
distances to raise the apparent physical inconsistency, is provided in Appendix A. 

In summary, we have the following possibilities for AE Cmut : 



• AE Cmut > for one possible £f™ h -value and AE Cmut < for the other. This is always 
true for intersecting footpoint segments and segments with "matching" footpoints. In 
this case we naturally select the £^ ch -value for which AE Cmut > 0. 

• AE Cmut > for both £^ cft -values. This can only be found in case of non-intersecting 
footpoint segments, where Cf^ h = J^f^f 1 - In this case, therefore, we naturally select 
this unique Cf^ h - value. 
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• AE Cmut < for both £^ c/i -values. This also happens exclusively in cases of non- 
intersecting flux-tube footpoint segments, and is because the connectivity-matrix cal- 
culation of Section 12.11 relies only on B z and is hence independent of a- values. One 
may envision an improved connectivity-matrix calculation in which the minimization 
functional of Equation (jSJ) is modified to include a-values - this will aim toward con- 
necting closely seated, opposite-polarity partitions with like-sign force-free parameters. 
At the present stage, however, we cannot physically accommodate negative free-energy 
increments, so we enforce Cf^ h = in these cases. Point taken, from practical expe- 
rience this assumption does not incur a large change in the free energy E Cmut as this 
energy term is dominantly influenced by flux tube pairs with intersecting or "matching" 
footpoint segments. 

Therefore, by means of Equation (123]) we both determine a unique £^ c/l -value for all 
off-diagonal helicity- matrix elements and ensure a symmetric helicity matrix. Hence, our 
pursuit for a minimum free magnetic energy has also led to a corresponding value for the 
relative magnetic helicity. 

Can we reach an even smaller, but still positive, value for the mutual term of the free 
magnetic energy E Cmut l The answer is yes, if we relax the assumption Cf^ se = 0. In case 
\^tm h \ > fc> r example, applying |£^ se | = 1 such that Cf^ h Cf^ se < 0, one may find 
smaller, and in some cases positive, increments AE Cmut . Moreover, even for cases where 
AE Cmut < for both C?^ h - values and Cf™ e = 0, we can always achieve AE Cmut > if we set 
£fm Se = 1 or Cf^ e = —1. All these possibilities are mathematically feasible. However, they 
lead to rather "exotic" physical situations of flux tubes winding around each other without 
necessarily intersecting footpoint segments. More importantly, they give rise to potentially 
uncontrollably high values of the mutual-helicity magnitude \H mmut \ or to helicity senses 
(chiralities) that run counter to the expected ones from observations of the active-region 
corona. We therefore follow the most conservative approach of keeping Cf^ se = and 
setting £?m h = in case it only yields a negative AE Cmut . In essence, our approach suggests 
that a realistic state for a non-potential, force-free magnetic configuration is achieved when 
the free magnetic energy is assumed minimum and the relative magnetic helicity is allowed 
to evolve in this respect, considering only arched and not braided flux tubes. Whether solar 
magnetic fields are in a minimum free-energy state is, of course, an open question. In Section 
[3] we show, however, that validating our approach with existing, generally accepted energy 
and helicity formulas leads to a fairly good agreement. 

Summarizing, given a collection of N discrete, slender, force-free flux tubes with flux 
contents $; and force-free parameters oti, where I = {1,2, ...N}, we write 
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(i) the total magnetic energy as E t = E p + E c , where E p is the potential magnetic energy, 
and E c is the free magnetic energy. The latter term is given by 



N N N 



E c = Ad 2 J2a^ x + ±J2 E (24) 



in 

1=1 1=1 m=l,ljim 



(ii) The relative magnetic helicity H m is given by 

N N N 



H m = 87Td 2 A + E E ■ ( 25 ) 

1=1 1=1 m=l,lj^m 

The scaling constants A, A are given by the least-squares best fit of Equation (JT5]) . the 
mutual-helicity terms Cf^ h are explained and calculated in Appendix A, and d is the size 
element. The force- free parameters on and flux contents of the flux tubes are inferred as 
described in Section 12.11 

It should be mentioned here that the methodology we describe is more general than 
Equations (1241) . (125]) in the sense that a nonzero Gauss linking number can always be ac- 
commodated in different physical settings. Moreover, our method can work for any given 
connectivity matrix. In this application, we physically favor both the connectivity matrix 
calculation of Section I2TT1 and C-f™ e = 0, thus reaching Equations (1241) . (12"5]) . 

A detailed error propagation analysis leading to the uncertainties 5E C and 5H m for E c 
and H m , respectively, is provided in Appendix B. 

In addition, one may infer the lowest possible free magnetic energy E CWT that corre- 
sponds to a given amount of relative helicity for the NLFF field. This is simply t he LFF 



1958; 


Tavlor 


1974. 


1986) 



we use the NLFF relative helicity H m inferred 
by Equation (|25|) and we calculate an effective constant a- value from Equation ( ITS"]) , i.e. 

Md 2 F lm E p ' 1 ' 

Substituting a from Equation (1261) into Equation (1T4"|) for the LFF free energy, then, we 
obtain 

tt2 tt2 
t? _ m rp _ m (97\ 

CWT (87rd) 2 F lin E p ' CWT (Md) 2 A^ ' 1 1 

due to Equation (1T5|) . The Woltjer- Taylor free magnetic energy E CwT of Equation (T27|) will 
be used as sanity check in the following; all calculated free magnetic energies E c must be 
larger than this limiting value. 
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2.3. A special case: potential-field configurations 

Our choice to calculate the magnetic connectivity matrix via simulated annealing (Sec- 
tion [2J]) implicitly assumes that the studied vector magnetograms are non-potential, that is, 
they include significant electric currents I and force-free parameters a for at least one of the 
major partitions identified. Although the non-potent iality of solar active regions is a long- 



known fact (e.g.. lZirin fc Wanglll993t iLeka et al.lll996l ). our methodology includes a physical 



inconsistency in case a potential-field configuration, observed or modeled, is subjected to 
the analysis: the simulated annealing method, favoring strong magnetic polarity inversion 
lines, provides a connectivity that is generally incompatible with that of the potential field. 
The question to ask, then, is how to determine whether a vector magnetogram is basically 
potential. 

One might argue that mutual-helicity terms £^ ch $/$ m algebraically cancel to zero in 
case of potential fields, along with individ ual self-helicity f actors T/ that tend to zero - 



otherwise, the relative-helicity expression of iDemoulin et al.l (120061 ) (Equation OH])) and our 



final formula (Equation (|25|) ) are not valid. In practice, however, an observed active- region 
magnetogram is not flux-balanced. This will inhibit an algebraic cancellation of mutual- 
helicity terms that may sum up to significant nonzero helicity values. At the same time, 
our expression for the free magnetic energy (Equation ( I24p ) will give values close to zero 
because a ~ for all partitions, and hence for all possible flux-tube connections. Hence, the 
physical inconsistency mentioned above leads to another inconsistency, namely to situations 
of near-zero magnetic free energy and strongly nonzero relative helicity. 

There are several methods to determine the degree of non-potentiality in observed active 
regions. We are currently in the process of identifying the most viable and practical of them, 
that will be reported in a future publication. One of these methods is to compare the observed 
horizontal-field components with those of the potential field obeying to the observed vertical- 
field distribution. Such a comparison may rely on scaling indices, correlation coefficients, 
and/or standard deviations between the observed and modeled components. Another method 
is to provide a flux-weighted mean a-value a and a respective uncertainty 5a from the a- 
values of all partitions. In case \a\ < n a 5a, where n a is a given significance level, the active 
region may be considered potential. Further, a flux-weighted mean magnetic-shear angle in 
the region is certainly another index of non-potentiality. In case of a nearly potential active- 
region magnetogram, rather than performing simulated annealing, one might perform a 
potential-field extrapolation and infer the connectivity matrix by line-tracing of the resulting 
magnetic field lines. Alternatively, one might stop the calculation at this point and set E c = 
and H m = for the active region of interest, thus saving computing time and avoiding large 
uncertainties owning to the magnetic-flux imbalance. 



-17- 



Method validation 



3.1. Benchmarking: volume formulas for magnetic energy and helicity 



Well-known formulas for magnetic energy and helicity can be used in case one knows the 
three-dimensional magnetic configuration above a two-dimensional planar boundary, where 
the magnetic field vector is fully known. For a three-dimensional flux-balanced magnetic 
structure the energy budgets are 



E, 



1 

8^ 



B 2 dV , E„ 



1 

8^ 



B 2 dV , and E c — E t — E p 



(2? 



where B — |B| is the magnetic field strength in the calculation volume V and B p 



IB, 



is 



the respective (also fully known) field strength of the potential field in V, where B and B p 
share the same normal component on the lower boundary. 



For the relative magnetic helicity H m in V. lBerger fc Fieldl (11984 ) and lFinn fc Antonsen 



(119851 ) derived two equivalent analytical forms valid for NLFF fields, namely 



H m = I (A±A p )-(B T B p )dV 
Jv 



(29) 



where A p and A are the vector potentials for B p and B, respectively. Although B p and 
B may be exactly known, proper knowledge of A p , and especially A, is a much more de- 
manding task. Point taken, the substantial value of Equation ( |29i) is that one may use two 
gauge-dependent, non-uniqu e expressions for A p and A to obtain a gauge-invariant, unique 



999). This p aved the way for implementatio n s suc h 



DeVore! (120001 ) and lLongcope fc Malanushenkci (120081 ) 



expression for H m (see also iBerger 
as the vector-potential expressions of 
among others. These expressions, however, are valid for the s emi- infinite space (half-space) 
above the lower-boundary plane. Recently, IValori et al.1 (120111 ) showed that if the vector po- 



tentials A p and A are corrected for the fact that the calculation volume V, wherein B p and 
B are known, is finite, then the resulting H m - values from Equation (1291) may be significantly 
different in amplitude and sometimes even in sign. In particular, these authors extended the 
analysis of [DeVoreJ (120001 ) and found that if V extends between (xi, x 2 ), (yi, 1/2), (21, z 2 ) in a 
cartesian coordinate system, then 



A r 



c — V x (z 



4>dz') 



(30) 



where is the scalar potential generating B p (B p = — V0) and c 



(c x ,Cy), with 



c x = -(1/2) flB Pz (x,y',z = z 1 )dy' 



(31) 
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and 



A = A p (x,y,z = zi) - z x / Bdz' , (32) 

J 21 

by choosing a gauge such that z ■ A p = z ■ A = everywhere in V. 



To test our NLFF expressions we will apply the derivations of IValori et all (120111 ) for 
A p and A to Equation (129|) for relative helicity in cases where we have performed NLFF field 
extrapolations on photospheric vector magnetograms. From these extrapolations and their 
potential-field counterparts we will also calculate the magnetic energy budgets of Equations 
(f2"8"|) . To be perfectly consistent with the purposes of this validation test, the connectivity 
matrix J^j on will not be inferred by simulated annealing (in this particular case only), but 
by tracing the NLFF-field-extrapolated lines. 



3.2. Validation results 



Validating our NLFF energy and helicity calculation method is a nontrivial exercise . 
At first glance, on e might rely o n analytical NLFF field models (e.g., lLow fc Lou! 11990 



Regnier et al.ll2005l ; lRegnierll201ll ). However, there are two shortcomings in this approach: 
first, in some of these models the vector potential is unknown, so one still needs to calculate 
gauge-dependent values of it to infer helicity. Second, and most importantly, the lower- 
boundary configuration is very simple in these models resulting in a very limited magnetic 
connectivity matrix. As we intend to emphasize the practical aspect of the methodology 
presented here, we choose to validate the method on real vector magnetograms, using gen- 
erally accepted energy and helicity formulas applicable to the coronal volume above the 
lower-boundary magnetogram (Section 13. ip . This volume is determi ned by a NLFF fiel d 
extrapolation. We use the cartesian optimization method developed by IWiegelmannl (120041 ). 
where the divergence of the magnetic field vector and the Lorentz force are simultaneously 
minimized for the configuratio n to converge to a NLF F state. No preprocessing of the 
boundary vector magnetogram (jWiegelmann et all 120061 ) was attempted in this case. 



Browsing through a sizable collection of vector magnetograms (jTziotziou et al.ll20121 ) we 
carefully selected 19 of them in which (i) the fractional magnetic flux imbalance is relatively 
low, thus allowing the majority of the unsigned flux to participate in the connectivity matrix, 
and (ii) the NLFF field extrapolation has worked well, with acceptable minimizations of 
the divergence of the field vector and the Lorentz force, and with a convergent solution 
exhibiting a total magnetic energy larger than the potential-field energy. To ensure that the 
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extrapolated three-dimensional field solution is valid we consider the dimensionless parameter 

IV • B\ 



V 



V3[(dB x /dx) 2 + {dBy/dy) 2 + (dBjdz) 2 ] 1 / 2 



(33) 



The flux- weighted mean of T> for all used magnetograms is < 0.1 indicating roughly divergence 
free, and hence valid, field solutions. The data were acquired by the IVM, with a binned 
pixel size of 2.2", and the Spectropola rimeter (SP; iLites et al.ll2008[ ) of Hinode's Solar Optical 
Telescope (SOT; iTsuneta et al.ll2008l ). SOT/SP is a spectrograph observing in two magnet- 
ically sensitive photospheric spectral lines; Fe I at 6301.5 A and 6302.5 A, respectively, with 
a spectral sampling of 21.6 mA . Full instrument resolution (~ 0.32") corresponds to a pixel 
size of ~ 0.158". The data used here, however, have been acquired in fast scanning mode, 
corresponding to ~ 0.308" per pixel. We have resolved the azimuthal 180° ambiguity in 
these dat a using the Non-P otential Field Calculation (NPFC) method of iGeorgoulisI (120051 ) 
- see also IGeorgoulisI (120121 ) . Then, to further expedite the extrapolations, we ha ve spatially 
binned the data to ~ 1.25" per pixel, indicating a spatial resolution of ~ 2.5". ILites et al. 



( 120081 ) reported line-of-sight and transverse-field uncertainties equal to 2.4 Mx cm~ 2 and 41 
Mx cm~ 2 , respectively, for the quiet Sun. In the following calculations we use uncertainties 
of 5 G and 50 G, respectively, for SOT/SP, and 50 G and 100 G, respectively, for IVM data. 
The same uncertainties apply for the calculations performed in the following Sections. 

Each vector magnetogram is then subjected to both volume- integral energy and helicity 
calculations (Section 13. 2j) and the NLFF surface calculations of this work. A crucial point 
here is that, for a direct comparison between the resulting energy and helicity budgets, one 
must use the connectivity matrix inferred by the NLFF field extrapolations. This is what 
we have done for this test only, inferring the various connectivity matrices by line-tracing 
the respective extrapolation results. The results of the comparison for the relative magnetic 
helicity and the free magnetic energy are given in Figures |3^i and [3b, respectively. Error bars 
correspond to uncertainties calculated by our NLFF field method while the red lines denote 
equality between the compared budgets. 

In terms of the relative magnetic helicity magnitude (Figure Et) we notice a fairly good 
agreement between the volume-calculated and the surface-calculated values - most points 
are within uncertainties from the equality line. The inferred helicity sense for extrapolations 
and our surface calculations agree in all 19 cases. For relatively small relative helicities 
(< io 42 Mx 2 ), our method appears to overestimate the relative helicity, albeit generally 
within applicable error bars. The linear (Pearson) and rank order (Spearman) correlation 
coefficients are significant, ranging between 0.63 and 0.74. 

The results of the comparison in magnetic free energy (Figure |3p) are similar to that of 
the relative helicity magnitudes. For small free energies (< 10 31 erg), our surface calculation 



-20 - 



seems to overestimate, within error bars, the respective values while the opposite, beyond 
error bars, occurs for free energies > 10 32 erg. The Pearson and Spearman correlation coef- 
ficients are also significant, higher than in the case of relative helicity, and ranging between 
0.74 and 0.75. We therefore conclude that the approximation of the free magnetic energy 
with the expression of Equation (j5J) is a reasonable choice, despite shortcomings, given the 
circumstances and the incomplete information. 

It is worth noting at this point that discrepancies between the volume and the surface 
calculation of magnetic free energy and relative helicity may be due to errors and ambiguities 
inherent to both our calculations, given the assumptions adopted, and the NLFF field ex- 
trapolations, including the calculation of the vector potentials that participate in the relative 
helicity formula of Equation (T2l?|) . Given the numerical methods involved, the derived vector 
potentials A p and A reproduce the magnetic field vectors B p and B, respectively, within 
non-negligible differences. 

In Figure [3b we have also plotted the Woltjer- Taylor free magnetic energy E CwT of 
Equation (I27p (crosses). These energies are to be viewed as sanity checks since the NLFF- 
field free magnetic energy E c cannot be smaller than them. In only one case in Figure [3b does 
E CwT exceed E c (fourth point from left). This, however, is within the applicable error bar. 
We further notice that some E CwT -va\ues are up to ~ 2 orders of magnitude smaller than 
E c , hence being unrealistically low, while some are quite close to their respective _E c -values. 
This may suggest the plausibility of the Taylor relaxation in at least some active regions. 
Although this discussion exceeds the scope of this work, we note in passing that the validity 
of the Taylor relaxation in the solar corona, according to which a relatively isolated magnetic 
configuration may relax in a state of minimum free (LFF-field) magnetic energy for a given 
magnetic helicity budget, is still a subject of debate. Perhaps calculation methods such as 
the one proposed here can provide further clues to judge the validity of this hypothesis. 

In brief, we conclude that our surface-based NLFF field calculation method manages 
to provide both magnetic free energies and relative helicities in fairly good agreement with 
generally accepted, but computationally much more intensive and model-dependent, volume- 
calculation techniques. It is, therefore, a viable method that can be exploited further and in 
larger data sets of solar vector magnetograms. 
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4. Results: NLFF field energy and helicity calculations 

4.1. NOAA ARs 8844 and 9165, in comparison with Paper I 

For the purposes of comparison with Paper I we use here the same two vector magne- 
togram timeseries, namely those of NOAA ARs 8844 and 9165, both acquired during daily 
observing cycles of the IVM on 2000 January 25 and September 15, respectively. The IVM 
data acquisition and selection and the properties of each AR are described in detail in Sec- 
tion 5.1 of Paper I. In brief, NOAA AR 8844 was a small emerging flux region, visible in 
the solar disk between 2000 January 24 and 27. No significant eruptive activity originated 
from this AR. On the other hand, NOAA AR 9165 was a complex, persistent, and eruptive 
region hosting a number of eruptive M-class flares (see Figure 2 of Paper I). Both ARs were 
fairly well flux-balanced within the IVM field of view: flux imbalance was < 5% for NOAA 
AR 8844 and < 10% for NOAA AR 9165 (Figure 3 of Paper I). The mean unsigned flux 
was ~ 5.1 x 10 21 Mx for NOAA AR 8844 and ~ 17 x 10 21 Mx for NOAA AR 9165, so with 
a unsigned-flux ratio ~ 3.4. Perhaps more relevant in this case is the total flux ^-f° n 
that participates in the connectivity matrices (Section 12. ip . where NLFF free energies and 
helicities stem from. The mean connected-flux values are ~ 3.6 x 10 21 Mx and ~ 12.8 x 10 21 
Mx for NOAA ARs 8844 and 9165, respectively, so with a connected-flux ratio of ~ 3.6. We 
notice that most of the unsigned flux for both ARs participates in the connectivity matrices. 
This is because both ARs are nearly flux-balanced. The little remaining flux that does not 
participate is either too disperse to be included in the flux partitioning or is judged to be 
connected with opposite-polarity flux concentrations seated beyond the field of view. 

Figure H] provides the timeseries of the calculated NLFF relative helicity H m (Figure H^i) 
and free energy E c (Figure Hb) for NOAA AR 8844. For an immediate comparison we have 
also plotted the respective LFF values from Paper I (gray curves and shades) while the NLFF 
field calculations of this work are shown with blue curves and shades. Assuming that within 
the observing interval of ~ 2 hours there was no significant change in H m and E c we define 
mean values H m and E c for H m and E c , respectively, accompanied by the respective standard 
deviations. We find H mNLFF = (1.18±0.45) x 10 42 Mx 2 and E CNLFF = (1.52±0.42) x 10 31 erg 
in the NLFF field approximation. Mean values are shown by the blue and gray straight lines 
for the NLFF and LFF calculation, respectively, while the respective standard deviations 
are shown by the blue- and gray-shaded areas. We notice that (i) the LFF and the NLFF 
approaches give values of H m and E c that are fairly close to each other (generally within 
uncertainties) at least for this small AR, and (ii) the NLFF field approximation gives an 
overall smoother evolution with smaller standard deviations. A summary of the values and 
uncertainties for the NLFF H m and E c is provided in Table [TJ 
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Figure [5] provides the timeseries of the calculated NLFF relative helicity H m (Figure 
[5^) and free energy E c (Figure [5b) for NOAA AR 9165. In this case the mean values 
H m and E c are higher than in NOAA AR 8844 (H mNLFF = (-7.3 ± 1.4) x 10 42 Mx 2 and 
E CNLFF = (5.3 ± 1.4) x 10 31 erg). The differences between the LFF- and the NLFF-field 
approximations in NOAA AR 9165 are also larger: in, general, the LFF-field approximation 
tends to overestimate both the magnitude of the relative magnetic helicity and the free 
magnetic energy. The overestimation caused by adopting the LFF-field approximation is 
~ 1.8 for H m and ~ 2.5 for E c . This is understandable and expected as the LFF-field 
approximation assigns a fixed force- free parameter a with a value determined by the strongest 
(most flux-massive) non-potential field configurations in the AR (see the analysis of deriving 
a single maximum-likelihood a- value in Paper I). In addition, the effect of obtaining smaller 
uncertainties for H m and E c in the NLFF-field approximation is more evident in NOAA AR 
9165. A summary of the mean values and respective uncertainties for NOAA AR 9165 is 
also provided in Table [TJ 

Comparing the mean values H m and E c for the two studied NOAA ARs 8844 and 9165 
we find a ratio of 6.2 ± 2.9 for H m and ~ 3.5 ± 1.3 for E c between the flaring and the 
non-flaring regions. These ratios are roughly similar to the unsigned- and connected-flux 
ratios (~ 3.4 and ~ 3.6, respectively) but smaller than those reported in Paper I for the 
LFF-field calculations of H m and E c . In that work, both ratios were ~ 9. We conclude 
that the overestimation of the magnetic free energy and relative magnetic helicity in the 
LFF field approximation is higher for larger, more complex active regions. Point taken, 
the flaring NOAA AR 9165 still shows much larger free-energy and relative-helicity budgets 
compared to the flare-quiet NOAA AR 8844. This quantitative disti nction between eruptive 



and non-eruptive active regions is studied by lTziotziou et al.l ( 120121 ). 



Concluding our calculations on NOAA ARs 8844 and 9165, we comment on the contribu- 
tions of the self and mutual terms to the free energy and relative helicity budgets (Equations 
f l24"|) . f l25|) . respectively) in the two studied ARs. We find that mutual terms overwhelmingly 
dominate these budgets: on average, for NOAA AR 8844 self terms contribute (0.1 ±0.07)% 
of the free magnetic energy and (0.1 ±0.2)% of the relative magnetic helicity. The respective 
percentages for NOAA AR 9165 are (0.5 ± .2)% and (0.4 ± 2)%. These findings are in 



qualitative agreement with previous works (IRegnier et al.l 120051 ; IRegnier &: Canfieldl 120061 ) 



and suggest that twist and writhe (contributing to self helicity) are numerically far less 
important than the mutual helicity caused by the interaction between different flux tubes. 
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4.2. NOAA AR 10930 



An appealing aspect of our analysis is that it can be applied to long-term timeseries 
of active-region vector magnetograms in order to reveal the temporal variation of magnetic 
energy and helicity budgets in the studied regions. The input vector magnetograms for this 
purpose should ideally exhibit constant quality and a fixed, high cadence. This is a big 
challenge for agin g ground-based mag netographs such as th e IVM. The Vector Spectromag- 
netograph (VSMi lHenney et al.ll2009l ) of the SOLIS facility (IKeller et al.ll2003l ) has achieved 
important advances on the quality front but, by design, it does not exhibit a cadence higher 
than a few hours. The space-based SOT/SP onboard Hinode exhibits unprecedented spatial 
resolution and a constant data quality due to the lack of atmospheric interference but, again, 
by design it only allows a cadence of a few hours. A lasting solution will have been achieved 
when the Helioseismic and Magnetic Imager (HMI; IScherrer et al.l 120 111 ) onboard the Solar 
Dynamics Observatory (SDO) releases constant- quality vec t or ma gnetograms of solar active 
regions with a fixed cadence of 12 minutes (see ISun et al.l (120121 ) for an example). At this 
time, however, for the purpose of showing the magnetic energy and helicity variations in an 
active region over a period of days, we present results obtained by processing a timeseries of 
SOT/SP vector magnetograms of NOAA AR 10930. 

NOAA AR 10930 appeared in the earthward solar hemisphere in 2006 December. It 
was an intensely flaring (and eruptive) region hosting ~ 50 C-class, 5 M-class and 4 X-class 
flares before rotating beyond the western solar limb. The AR has been studied in extreme 
detail by dozens of works; in some of them estimates of the magnetic free energy and the 



relativ e magnetic helicity in the region have been published. For example, iRavindra et al. 



( 1201 ll ) reported that, by 2006 December 13, when the AR hosted a X3.4 flare, more than 
—6 x 10 43 Mx 2 of relative h elicity had b een in jected in the AR. A slightly more conservative 
estimate was published by iPark et al.l (120101 ). with a relative helicity ~ —4.3 x 10 43 Mx 2 



before the flare. As to the free ma gnetic energy of the AR, iHe et al.l ( 120111 ) estimated it 



within (1.25 — 1.4) x 10 33 erg, while iGuo et al.l (120081 ) reported a value of ~ 2 x 10 32 erg. 



We have selected 30 vector magnetograms of the region acquired by Hinode's SOT/SP 
between 2006 December 8 and December 14. The heliographic vertical magnetic field compo- 
nent of six of them is shown in Figure [6] with Figures [6^ and[6j: corresponding to the first and 
last magnetogram of the series, respectively. The magnetograms were first disambiguated 
using the NPFC method and then were spatially binned to ~ 0.62" per pixel, or to a spatial 
resolution of ~ 1.23". 

The evolution of the relative magnetic helicity in NOAA AR 10930 is shown in Figure [T^l. 
A distinct feature of the pre-emption evolution in the AR is that, before 12/10, the relative- 
helicity budget is rather low, of the order (— 5± 3) x 10 42 Mx 2 . Over the next two days (12/11 



-24- 



- 12/12), however, the helicity budget increases drastically to reach ~ — 1.1 x 10 43 Mx 2 . This 
peak coincides with a cluster of C-class flares (see the respective GOES 1-8 A X-ray flux 
in Figure [7U) that appear to be eruptive, as can be judged by the repetitive Type-II bursts 
recorded in the frequency-time radio spectra of the WAVES instrument onboard the WIND 



mission (B ougeret et al. 



Nelson fc Melrose 



1995). Type II activity implies shock-fronted CME occurrences 
1985h . Although neither GOES nor WIND/WAVES have spatial 



(e.g., 

resolution, there is little doubt that the observed eruptions originate from NOAA AR 10930, 
as it was the only AR present in the visible solar disk at the tim e. The CMEs and their 
locations are also confirmed by the SoHO/LASCO CME database (jYashiro et al.l 120041 ) . 



Around the start of 12/12, perhaps due to eruptive activity, the helicity budget appears 
generally smaller, of the order — (7 ± 3) x 10 42 Mx 2 . Hours before the X3.4 flare, early on 
12/13, however, the relative helicity budget increases substantially to exceed ~ —1.5 x 10 43 
Mx 2 . At this time of peak helicity the eruptive flare occurs, accompanied by a fast halo CME. 
Immediately after the eruption, late on 12/13, the AR appears to have lost ~ —5 x 10 42 Mx 2 
of helicity, perhaps in the CME. For a period of ~ 12 hours until the end of the observing 
interval the helicity budget appears to be of the order ~ — 1 x 10 43 Mx 2 . 

Around the time of the X-class flare we calculate a r elativ e magnetic helicity budget that 
is a factor of ~ 3 lower than the estim ate of iPark et al.l (120101 ) and a factor of ~ 4 lower than 
the estimate of iRavindra et al.l ( 1201 ll ). Our lower helicity estimate is consistent with (i) the 
fact that we calculate helicity from closed-field connections only, thus using only a fraction of 
the unsigned flux (a small fraction in this the AR shows significant flux imbalance - 

Figure [Tt) and (ii) our methodology, that minimizes the free magnetic energy first, and then 
keeps a consisten t helicity magnitu de. In addition, our study is qualitatively consistent with 
the assessment of IPark et al.l ( 120101 ) that helicity magnitude abruptly increases on 12/10 and 
thereafter. A detailed discussion of the process exceeds the scope of this work that mainly 
focuses on the proposed calculation methodology. 

The evolution of the magnetic energy budgets in NOAA AR 10930 is shown in Figure 
[7b. Here we notice that while the potential energy is linearly increasing over the observing 
interval, in agreement with the increase of the unsigned flux (Figure Et), the free magnetic 
energy increases with a faster rate on 12/10 and thereafter, when the relative helicity magni- 
tude increases. The peak helicity at the time of the multiple eruptive C-class flares coincides 
with a free energy of ~ 2 x 10 32 erg. The free energy is kept at approximately these levels, 
albeit showing a moderate increasing trend, until late on 12/12. Thereafter, it increases 
significantly to peak at ~ 5 x 10 32 erg around the time of the X-class flare. After the flare, 
following the decrease of the helicity budget, the free energy decreases to ~ 3 x 10 32 erg, to 
remain at these levels until the end of the observing interval. Given the small uncertainties, 
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the decrease of ~ 2 x 10 32 erg in the course of the eruption appears significant. 



Our free-energy estimates agree qualitatively w ith those of iGuo et al.l ( 120081 ) but they 
are significantly lower than those of iHe et al.l ( 120 111 ). This is in agreement with our effort 
to keep the free energy minimum, at the same time excluding from calculation all magnetic 
connections that close beyond th e field of vi ew. T his being said, the postflare free energy 
decrease is quite consistent with iHudsonf s (120111 ) assessment that X-class flares typically 
dissipate ~ 10 32 erg of magnetic energy. 

In summary, we find that a plausible physical interpretation of the dynamical evolu- 
tion of the eruptive NOAA AR 10930 can rely on our calculation of the relative magnetic 
helicity and free magnetic energy in the region. Moreover, the estimated magnetic energy 
and relative helicity budgets are consistent with the lowest published estimates, as expected. 
Given that our calculations (i) stem from a unique solution for the magnetic connectivity ma- 
trix, (ii) do not depend on magnetic field extrapolations, and (iii) are relatively inexpensive 
computationally, we argue that the proposed method is both viable and practical. 



5. Discussion and Conclusions 



Motivated by the need to achieve a practical, realistic, and self-consistent assessment of 
the magnetic energy and relative magnetic helicity budgets in solar active regions we first 
tackled the problem in the simplified LFF field approximation (IGeorgoulis fc LaBontdl2007 
- Paper I). In this case the solution was unique and dependent on the single, fixed force- free 
para meter used. However, it is known tha t there are multiple ways to infer this single a-value 



U.e. jLeka fc Skumanichlll999t lLekalll999l ) and the one used in Paper I was but one of several 
methods. Different a-values give different solutions for the magnetic energy and helicity 
budgets. Moreover, the LFF fie l d app roximation is generally unrealistic in active-region 
scales (see, however, iMoon et al. (booA for a different view). For this reason the analysis 
of Paper I was conceived as the first step toward a more realistic approach of performing 
magnetic energy and helicity calculations. This step, relying on a NLFF field approximation, 
is taken in this work. 

Multiple methods to calculate the magnetic energy and relative helicity budgets of solar 
active regions assuming NLFF magnetic fields are long present in the literature. Virtually 
all of them rely either on a photospheric velocity flow field o r on a three-dimensional NLFF 
field of the active-region corona (for a review, see iDemoulinl (120071 ) and references therein). 
Both the photospheric flow field and the three-dimensional coronal NLFF field, however, are 
not uniquely defined. Therefore, the resulting NLFF-field energy and helicity budgets are, 
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again, model-dependent. An apparently more robust method to cal culate the spinning and 



braiding helicity in observed solar active regions was intr oduced by lLongcope et al.l (120071 ) 



and applied to NOAA AR 10930 by lRavindra et al.l ( 1201 ll ). The method also uses a velocity 



field but now this field is obtained by tracking the motions of photospheric flux partitions. 
These partitions are inferred as described in Section 12.11 and Figure [TJ Feature tracking on 
individual partitions should be more stable than calculating the entire flow field, although 
limitations have been reported for this well. 

Here we follow a different approach depending on neither photospheric flow fields nor 
the unknown coronal three-dimensional field. Our method depends on a lower-boundary 
(photospheric or chromospheric) magnetic connectivity matrix that can be inferred either 
by a NLFF field extrapolation or otherwise. Had we used an extrapolation, our results would 
also be model-dependent. Although our method is general enough to accommodate any con- 
nectivity matrix, we propose and use a unique connectivity-matrix solution for a given flux- 
partition map. This solution relies on a simulated-annealing algorithm designed to minimize 
the distances of connected opposite-polarity partitions, thus emphasizing strong polarity in- 
version lines. This connectivity methodology has been s uccessful in distinguishing flaring 



from non-flaring active regions (IGeorgoulis fc Rustl 120071 ) and will be shown to be of fur 



ther such importance in a much larger statistical sample (work in preparation). Besides the 
boundary connectivity, our method alleviates the need for flow fields and three-dimensional 
coronal field vectors by, first, minimizing the magnetic free energy of the active-region corona 
and, second, keeping the dominant mutual-helicity terms consistent with the global (within 
active-region scales) energy-minimization principle. Self terms of magnetic energy and he- 
licity are calculated by generalizing the LFF field analysis of Paper I, whi le mutual terms 



are ca lculated by a practical implementation of the method introduced by iDemoulin et al. 



(120061 ). 



To validate the proposed method we use connectivity matrices derived from NLFF field 
extrapolations because comparison is then based on the NLFF model-dependent energies and 
helicities. For the validation we use real active-region magnetograms, thus avoiding synthetic 
NLFF fields with a simpler lower boundary and hence a smaller and cruder connectivity 
matrix. At the same time, the approximate validity of the volume-calculated energies and 
helicities is guaranteed by the use of well-known and accepted energy and helicity formulas. 
We find (Section 13. 2p that the results of known volume formulas that require a detailed 
three-dimensional field configuration are generally reproduced by our surface formulas that 
use only the connectivity matrices inferred from the extrapolations. This justifies the use of 
the free-magnetic-energy formula of Equation (l24p as a lower limit of the true free energy, 
in spite of its weakness to fully describe space-filling, force-free fields (Section 12 . 2 j) . 
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Following validation, we apply our method to the same set of active regions with Paper 
I and compare the results (Section 14.11) . We find them to be consistent, in general, but 
with the LFF field approximations overestimating the free magnetic energy and relative 
magnetic helicity budgets in such a way that overestimation is higher for larger, more complex 
active regions. In the NLFF field approximation the ratios of the free magnetic energy and 
the relative magnetic helicity between the two ARs are roughly similar to the ratio of the 
unsigned flux participating in the magnetic connectivity matrices (J&mtx ~ see Table [TJ, 
as opposed to ~ $ 2 , i.e., the square of the unsigned magnetic flux, obtained by the LFF 
field approximation of Paper I. In addition, we find that mutual-energy and helicity terms 
overwhelmingly dominate the respective budgets in both active regions with contributions 
in excess of 99.5%. 

Our NLFF field approach is then applied to a timeseries of vector magnetograms of the 
eruptive NOAA AR 10930, observed over a perio d of ~ 6 days (Section 14. 2p . The results 
corroborate the findings of iRavindra et al.l (120111 ) that the relative helicity in the region 
increased abruptly within ~ 24 hours, resulting in significant left-handed helicity in the 
region. The peak hel icity magnitude we find, however, is ~ 4 time s smaller than that of 



Ravindra et al.l ( 120111 ) and ~ 3 times smaller than that of iPark et al.l (J2010J). We also find 
that (i) an initial helicity magnitude decrease is associated in time with a cluster of eruptive 
C-class flares, and (ii) a more abrupt helicity decrease of ~ —5 x 10 42 Mx 2 relates closely in 
time with the eruptive X3.4 flare that climaxes the eruptive activity in the region over the 
observing interval. A similar decrease of ~ 2 x 10 32 erg in the region's free magnetic energy 
was also calculated at that time. Both the relative helicity and the free energy remained at 
these lower-budget levels until the end of the observing interval. To firmly establish these 
findings, however, more vector magnetogram data are necessary. 

The practical energy and helicity calculations in observed solar active regions being 
the scope of this work, the ultimate objective, also posed in Paper I, is to acknowledge 
and outline the possible role of mag netic helicity in the tr iggering of solar eruptions. We 
take a first step in this direction in iTziotziou et al.l (120121 ). where the method introduced 
here is applied to 162 active-region vector magnetograms to yield the first energy-helicity 
diagram of solar active regions. This diagram demonstrates a monotonic correlation between 
the free magnetic energy and the relative magnetic helicity in active regions, at the same 
time showing a segregation between flaring/eruptive regions and non-e ruptive ones. This 



finding reinforces the resul t s of previous works in o bserved active re gions (jNindos fc Andrews 



2004; 


LaBonte et a 




2007 


- see also 


Nindos ( 


2009) 



20071 ; iGeorgoulis et al.ll2009f ) and in theory (jZhang fc Low 



2001 



2003 



Nindos ( 20091 ) for a review) that eruptions leading to CMEs effectively transfer 



excess magnetic helicity from the Sun outward to the heliosphere. One is tempted to assert, 
therefore, that if CMEs are means to relieve the Sun from its excess helicity, then helicity itself 
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may play a key role in solar eruptions. Other than the heli cal kink instability, kn own to lead 
to eruptions in some observed filament de stabilizations (IRust k, LaBontd 120051 ). however, 
and the helicity- annihilation mechanisms of Kusano et al. (12003) that remain s to be proved, 
the alleged role of helicity is unknown. Only very recently, iKliem et al.l (120121 ) demonstrated 
that weakly kink-u nstable magnetic co nfigurations can represent observed solar-eruption 
features. Moreover, iRaouafi et al.l (120101) presented evide nce that several X-ray jets observed 
by Hinode's X-Ray Telescope (XRT; iGolub et al.l 120071 ) are preceded my micro- sigmoids, 
thus elucidating the possible role of the helical k ink instability even in small-scale eruptive 
activity. Finally, iPatsourakos fc Vourlidasl ( 120091 ) stereoscopically observed a small transient 
sigmoid that erupted giving rise to an observed EIT wave. In view of the above and other 
results, a detailed investigation of the role of helicity in the pre-eruption configuration of 
eruption-prolific s olar active regions is well justified. Such an investigation was sketched by 



Georgoulisl ( 1201 ll ) and will be the subject of a study currently in preparation. 
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A. Calculation of possible C™ c ch - values for two discrete magnetic flux tubes a 

and c 



Per iDemoulin et al.l (120061 ) . the mutual-helicity parameter C ac = C^° se + C^ h of a 
pair of flux tubes a and c can be calculated by progressively bringing the interacting pair 
from infinity to its prescribed position and geometry. This can be generalized for any set 
of discrete flux tubes, where (a, c) is now any given pair of tubes belonging to the set. For 
each case of interaction there are only two possible values of C^ r c ch that depend on (i) the 
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geometry of the pair, as reflected by the line segments formed by the footpoint locations of 
each tube, and (ii) whether tube a is "above" tube c (C^ r c ch = C- a ^ h ) or tube c is above tube 
a (J~-ac Ch = ^a6 )■ By stating that a tube is "above" another tube we imply that its apex is 
higher than th e appex of its m a te, wi th respect to the anchoring boundary that defines the 
open volume. iDemoulin et al.l ( 120061 ) further demonstrated a practical way of calculating 
C^ r c ch , namely by means of interior angles of the triangles formed by footpoint segments on 
the anchoring boundary plane. In the follo wing we calculat e L a J n ch for all cases pertinent to 
our study, both reproducing the values of IDemoulin et al.l (120061 ). and deriving values for 
cases that were not examined by these authors. 

In practice, C^ r c ch is the mean angle by which each line segment "sees" the other, normal- 
ized over 7r. Therefore, assuming a flux tube a with positive- and negative-polarity footpoints 
a + and a_, respectively (so its footpoint segment is a + a_), each footpoint "sees" the dipole 
c (with footpoint segment c + c_) by angles a a+ and a a _, respectively. Then, in case a is 
"above" c we have 

Q r c ch = + a a _) . (Al) 

Similarly, footpoints c + and c_ of tube c will "see" the segment a + a_ with angles a c+ and 
a c _ , respectively. In case c is "above" a we have 

£Z ch = ^(a c+ + a c _) . (A2) 

It should be mentioned that interior angles a x± (x = {a, c}) do not necessarily follow the 
trigonometric (counterclockwise) convention, so they can be positive or negative. Moreover, 
it is clear from Equations flH) and (|A2]) that |£^ c/l | < 1. 

All possible footpoint-segment configurations in our calculations appear in Figure 
Cases o f non-intersecting (Figu re Ek) and intersecting (Figure Eb) segments were also exam- 
ined by IDemoulin et al. J2006h . Given that our connectivity matrix has been constructed 
by connecting magnetic partitions, however, it is very common to find multiple connections 
connecting a given partition with others. Since all connections are viewed as slender flux 
tubes originating from the partition's flux- weighted centroid (Figures QJ>) we have cases 
of "matching" footpoints, as well (Figures [Sb, EH)- The apparent conflict with the principle 
that magnetic field lines cannot intersect may be raised by clarifying that these "matching" 
footpoints are, in fact, distinct footpoints but with distances falling into unresolved length 
scales within a given partition. Therefore, for the sake of simplicity they are thought to orig- 
inate from the same well-known location, i.e., the partition's flux- weighted centroid. This 
introduces some modifications in £^ cft -values and we calculate these modifications here. 

The four possible locations of footpoint-segment geometry, as illustrated in Figure El 
are as follows: 
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CASE A: Non- intersecting footpoints (Figure E^). It can then be found that the footpoints of 
flux tube a "see" flux tube c by the angles 



a 



a+ 



9, 



e, 



c + a + 



and a a _ = 8 c+a _ - 9 c _ a 



where 9 x±y± (x,y = {a,c}) are the azimuth angles of the segments y±x± (that is, with 
the trigonometric-circle origin at x±). Correspondingly, the footpoints of flux tube c 
"see" flux tube a by the angles 



Or 



-0 r 



>, : ' ' i ' (/_(_ "I - 9 c+a _ and ol c _ 9 ca 
Obviously, then, from equations ( lAlj) . (1A2j) we obtain 

parch parch 

^ac ^ac i 



e, 



(A3) 



reproducing iDemoulin et al.l ( 120061 ). In essence, in case of non- intersecting footpoint 
segments the two possible £^ c/l -values collapse to a single value, regardless of appex 
heights for the two flux tubes. 

CASE B: Intersecting footpoints (Figure [8b). In case a is "above" c (upper configuration) we 
find 



a 



a+ 



-e, 



c+a+ 



+ 9, 



c—a+ 



and a n 



2tt - 6 r 



0, 



In case c is "above" a (lower configuration), we further have 

a c+ = -9 c+a+ + 9 c+a _ and a c _ = -9 c _ a _ + 9, 
From Equations (1A1I) and (1A2[) . then, we find 

■•arch parch i 



c_a + ; 



part 
l —ac 



pare 
l —ac 



(A4) 



also reproducing IDemoulin et al.l (120061 ). 
CASE C: Positive "matching" footpoints (Figure [Sb). In this case we obviously have 

OC a , 



To calculate the remaining interior angles a a _ and a c _ we first notice that the formed 
triangle dictates 



+ a„_ + \a c 



71 



(A5) 



The orientation (sign) of a a _ , a c _ depends on the orientation of the angle (p ac between 
segments a and c. Angle a a _ always follows the sign of (p ac while angle a c _ always 
shows opposite orientation. Depending on the orientation of the triangle as a rigid 
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shape, one may further find from trigonometric analysis that the magnitudes of ct a _ , 
a c _ are given by 

|a _| = \9 c _ a _ - 9 a _ a+ ± ti\ 

c _ a+ -e c _ a _\ (A6) 



C -a + - Oc-a- =F 27r| 



It then becomes trivial to find which of the six possible cases applies for |a _|, \ot c _\ 
so that Equations ( 1A6j) satisfy Equation (IA5j) . The orientation of a a _, a c _ is then 
found when the respective orientation of <j) ac is found, that is, by determining whether 
one forms an interior (hence < 7i~) angle from a to c following (+) or opposing (-) the 
trigonometric convention. For the example of Figure [St we obviously have ac > 0, so 
a a _ > and a c _ < 0. 

Calculating the £^ c/l -values in this case we always find 

CZ ck C arc h < j (A7) 

which makes it straightforward to determine the preferred value so that AE Cmut > in 
Equation fl23|) . 



CASE D: Negative "matching" footpoints (Figure [HJd) . In this case we have 

«a_ = ol c _ = . 

The analysis is symmetric to that of Case C in the sense that a a+ always follows the 
orientation of <j) ac and a c+ always opposes it. From the formed triangle we have 

\<p ac \ + \a a+ \ + \a c+ \ = 7T , (A8) 
and a trigonometric analysis for different triangle orientations shows that 

\9 c+ a„ ~0 c+a+ \ (A9) 
|^c+a_ — Q C+ a+ =F 27r| . 



1 1 




1^. 






{ 



We solve the systems of Equations (1A8[) . (1A9j) to determine |a a+ |, |a c+ | and we infer 
the orientation of 4> ac to determine their orientations. In this case, as well, we find 

£Z ch £:? h < , (A10) 

so determining the desired value for AE Cmut > in Equation (1231) is straightforward. 



-32 - 



B. Uncertainties for Magnetic Free Energy and Relative Magnetic Helicity 

Assuming that the uncertainty in the calculation of the potential magnetic energy E p 
is negligible, the uncertainty 5E t of the total magnetic energy equals the uncertainty 5E C of 
the free magnetic energy E c . Hence 

5E c = 5E t . (Bl) 

As the free magnetic energy comprises of self and mutual terms, E Cse[f and E Cjnut , respectively, 
the uncertainty SE C depends on the uncertainties of these two terms, i.e. 

5E C = sJ{8E Cself f + (SE^J* , (B2) 
assuming, of course, that E Cself and E Cmut are independent of each other. 

Likewise, the uncertainty SH m of the relative magnetic helicity depends on the respective 
uncertainties of the self and mutual helicity terms, i.e. 

5H m = ^/(5H mself y + (5H mmut )i . (B3) 



B.l. Uncertainties in the self terms of free energy and relative helicity 

From Equation ( |T6|) . the uncertainty SE C(a) in case of a single force- free flux tube (that 
is, without mutual terms) depends on the uncertainties 5a, 5A, and 5X of the force-free 
parameter a of the tube and the fitting parameters A and A, respectively, of Equation (|T5|) . 
assuming that the tube's flux content $ is known without uncertainty. By standard error 
propagation and assuming that a, A, and A are independent of each other, we find 

^ C(s) =E C(s) [4(^) 2 + (^) 2 + 4(ln$) 2 (5A) 2 ] 1/2 ■ (B4) 

From Equation (1B4|) and assuming a collection of N flux tubes, the uncertainty SE Caelf in 
the self term of the free magnetic energy (Equation ( 1TTI) ) will then be 



\ 



N 



E^X/) 2 ■ ( B5 ) 
/=1 



Likewise, for a single force-free flux tube the relative magnetic helicity is given by 
Equation (fT9|) . Therefore, the uncertainty 8H m , s) in this case is 

5H m(s) = |^ m(s) |[(^) 2 + (^) 2 + 4(ln$) 2 (5A) 2 ]^ . (B6) 
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Then, from Equation (120]) in case of a collection of N flux tubes, the uncertainty 5H mself in 
the self term of the relative magnetic helicity is 

N 

In both Equations (1B5I) and (1B7I) each uncertainty 5a / of the involved flux tubes' force- 
free parameters is calculated by Equation (j3J), thus involving the uncertainties of the flux- 
weighted a-values of the connected partitions. 



5H„ 



'self 



\ 



B.2. Uncertainties in the mutual terms of free energy and relative helicity 

From Equation §ZQ, the mutual term of the free energy for a given pair (l,m) of flux 
tubes is given by 

E CmuU im = ^-a lm £Z h $i®m ■ (B8) 
The uncertainty SE Cmutt i m is then given by 



,5aim, 2 , ,5C arch 

01 lm *-lm 



$ E c mut ,lm = E CmutM i /( -) 2 + ( ra l T c h ) 2 ' ( B9 ) 



ignoring dependencies between ai m and Cf m ch for simplicity. Now assuming that each mutual- 
energy term E Cmut j m is independent from the other terms, the overall uncertainty SE Cjnut in 
the mutual term of the free magnetic energy is 



6E r 



N N 



\ E E i^muUlm? , (BIO) 
\ 1=1 m=l, Z^m 



for a collection of N flux tubes. Therefore, the problem becomes equivalent to determining 
the uncertainties 5ai m , where ai m = (l/2)(azi+a m ), and 5£? m ch , for each mutual-energy term 
(Equation (jB9|l ). 



For 6 aim we trivially have 



ly/(5 ai ) 2 + (5a m )* , (BH) 



5 ai m 

assuming that ai and a m are independent of each other. The uncertainty SC^ m ch of C-^ m ch 
depends sensitively on the uncertainty Sai m - In particular: 
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In case one of the two possible AE^^-values in Equation ff23l) is positive but the sign 
of ai m is ambiguous (5ai m > \ai m \), any of the two possible Cf m ch - values could be used. 
In case the sign of ai m is known with certainty (5ai m < \cxi m \) then the choice of one 
of the two possible £^ c?t -values is definitive. Therefore, 



SC arck I \Cf^ h -q m ch \ ; if AE<£ ut AE?L < and 6a lm >\a lm \ 
'"' _ " S " if AEj^AEgL < and (Ja, ro < |a Jfl 1 " 



In case both possible AE Cmut -vahies in Equation (1231) are either positive or negative 
but the sign of ai m is ambiguous (5ai m > |atj m |), then either the selected £? m ch - value or 
£^ c/l = could be used. In case <5«z m < |«; m |, the choice of the non-zero or zero £f m ch 
is definitive. Because this situation can only occur in case footpoint segments are not 
intersecting, so £^ c/i = £f m ch , we define an uncertainty 



■•arch ) l^lm 



,^£„r,-/, j \£7m h \ ! *7 AE^AE^t > and <$a lfn > |a Im| 



When the mutual term of the relative helicity is concerned, for a given pair (/, m) of 
flux tubes we have 

H mmKtM = CZ h ^m ■ (B14) 

Assuming, as usual, that <& m are known without uncertainties, the uncertainty 5H mmut j m 
is given by 

fiHnimutM = \H mmut ,im\5C? m ch ■ 
Then, the overall uncertainty 5H mmut j m is trivially given by 



6 Urn 



N N 



V E E ( tffl »W.«n0 2 > ( B1 6) 
^ 1=1 m=l,lj^m 



where each uncertainty SH mmut i m is calculated by Equation (IB 1 5 [) . 
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Fig. 1. — Example magnetic connectivity in NOAA AR 10254, observed by the IVM on 2003 
January 13, at 21:44 UT. Images (a) and (b) show the vertical magnetic field component 
in grayscale, saturated at ±1000 G with the contours bounding the identified magnetic 
partitions. Image (c) shows the respective flux-weighted a- value for each partition. The flux 
tubes identified by the magnetic connectivity matrices are represented by line segments of 
different colors and thicknesses, connecting the flux-weighted centroids of the respective pair 
of partitions. Only closed connections within the field of view are shown and considered; 
disconnected partitions are exclusively linked to opposite fluxes beyond the field of view. Red, 
cyan, and yellow segments denote magnetic flux contents within the ranges [5 x 10 19 , 5 x 
10 20 ] Mi, [5x 10 20 , 10 21 ] Mx, and [10 21 , 5x 10 21 ] Mx, respectively. The connectivity matrices 
shown enclose (a) ~ 10 22 Mx for the simulated-annealing and (b) ~ 8.3 x 10 21 Mx for the 
potential-field connectivity. Tic mark separation is 20". North is up; west is to the right. 
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Fig. 2. — Relation between the "scaled" potential energy FunEp and the square $ 2 of 
the unsigned magnetic flux for 56686 active-region line-of-sight magnetograms recorded by 
SoHO/MDI between 1996 and 2005. All these regions are located within ±30° from the 
solar central meridian at the time of observation. The red line is the least-squares best fit 
described by Equation (TT5T) . 
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Fig. 3. — Comparison between our NLFF surface-calculated values (ordinate) and well- 
known volume-calculated values (abscissa) for (a) the relative magnetic helicity (|-ff m |) and 
(b) the free magnetic energy (E c ) budgets. Equality between the two budgets is denoted 
by the red lines. In (b), crosses denote the Woltjer- Taylor free magnetic energies E CwT 
(Equation (I2"T)) ) linked to the respective free magnetic energies E c by blue dotted lines. 
Error bars (only to higher values in some cases, as including the lower-value error bars 
would result in negatives in these cases) correspond to our NLFF surface-calculated values 
for both plots. 




Fig. 4. — Comparison between the LFF and NLFF field approximations in calculating (a) 
the relative magnetic helicity budget H m and (b) the free magnetic energy budget E c in 
NOAA AR 8844, observed by the IVM for nearly two hours on 2000 January 25. The mean 
values H m and E c are indicated by straight lines (blue for the NLFF-, gray for the LFF- 
field calculations) while the respective standard deviations are indicated by the blue- and 
gray-shaded areas. 




Fig. 5. — Same as in Figure H] but for NOAA AR 9165, observed by the IVM for nearly five 
hours on 2000 September 15. 
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Fig. 6. — Heliographic vertical field components of selected vector SOT/SP magnetograms 
of NOAA AR 10930 binned to ~ 0.62" per pixel. The dates and universal times of the 
start of observations for each magnetogram are also shown. The vertical field component is 
saturated at ±2000 G.Tic mark separation is 20". North is up; west is to the right. 
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Fig. 7. — Temporal evolution of NOAA AR 10930, observed over a period of a few days in 
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2006 December by Hinode's SOT/SP, as reflected on (a) the calculated relative magnetic 
helicity budget, (b) the calculated magnetic energy budgets (with the total, potential-field, 
and free, magnetic energy shown by red, green, and blue curves, respectively), and (c) 
the unsigned magnetic flux (red) and flux imbalance (purple), including the unsigned flux 
participating in the connectivity matrix (blue). Also shown for reference are the respective 
(d) GOES 1-8 A solar X-ray flux, showing flaring activity, and (e) WIND/WAVES frequency- 
time radio spectrum, with Type-II activity indicating shock-fronted CMEs. All flares and 
CMEs originate from the region because it was the only AR present in the solar disk at the 
time of interest. 
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Fig. 8. — Possible geometrical positions of the photospheric footpoints of two discrete, arched 
flux tubes a and c: (a) non-intersecting footpoint segments, (b) intersecting footpoint seg- 
ments, including cases where tube a is "above" tube c (upper sketch) and where tube c is 
"above" tube a (lower sketch), (c) "matching" positive-polarity footpoints, and (d) "match- 
ing" negative-polarity footpoints. In all sketches, footpoint segments are depicted by thick 
black lines, while all connecting lines are cyan. All in terior angles of the for med triangles 
are also shown. Cases (a) and (b) were first studied by iDemoulin et al.l ( 120061 ) but they are 
reproduced here for reasons of completeness. 
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NOAA AR 


$ 


&MTX 


E p 


E c 


H m 




(xlO 21 Mx) 


(xlO 21 Mx) 


(x 10 32 erg) 


(xlO 32 erg) 


(xlO 42 Mx 2 ) 


8844 


5.1 ±0.2 


3.6 ±0.2 


2.97 ±0.01 


0.152 ±0.04 


1.18 ±0.5 


9165 


17.1 ±0.8 


12.8 ±0.4 


9.52 ±0.5 


0.53 ±0.14 


-7.3 ± 1.4 


Ratio 


3.4 ±0.2 


3.6 ±0.2 


3.2 ±0.2 


3.5 ± 1.3 


6.2 ±2.9 



Table 1: Synopsis of the mean unsigned magnetic flux ($), connected flux (i.e. unsigned 
flux participating in the magnetic connectivity matrix; <&mtx), current-free magnetic energy 
(E p ), free magnetic energy (E c ), and relative magnetic helicity (H m ) budgets for NOAA 
ARs 8844 and 9165 as obtained by our NLFF methodology. The third row refers to the 
ratio I-P9165/-P8844I between a given parameter P 9 i 65 of NOAA AR 9165 and the respective 
parameter Pg844 of NOAA AR 8844. In terms of mean free magnetic energy E C) ~ 99.9% 
for NOAA AR 8844 and ~ 99.5% for NOAA AR 9165 are contributed by the mutual terms. 
For the mean relative magnetic helicity H m the mutual-term contributions are ~ 99.9% for 
NOAA AR 8844 and ~ 99.6% for NOAA AR 9165. 



